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ABSTRACT 

It is easy to construct classical 2-state systems illustrating the behavior of 
the short-lived and long-lived neutral K mesons in the limit of CP conser- 
vation. The emulation of CP violation is more tricky, but is provided by 
the two-dimensional motion of a Foucault pendulum. Analogies are drawn 
between the pendulum and observables in neutral kaon decays. An emu- 
lation of CP- and CPT-violation using electric circuits is also discussed. 

I. INTRODUCTION 

Two-state systems abound in quantum mechanics, and have the virtue of permit- 
ting exact solutions [l]. One such system is that of the neutral kaons K° and K° 
||, which mix with one another || to form short-lived and long-lived mass eigen- 
states. Previously [|J we have shown that one can emulate this system using coupled 
resonant circuits, with the violation of CP invariance || exhibited in kaon decays 
corresponding to asymmetric coupling between the circuits. 

We sought to implement the suggestion of Ref. 0] in a physical device, such as 
a pair of coupled pendula which may be used to illustrate the CP-conserving limit 
||. In the course of this activity, we came upon a familiar system which has many of 
the features of the neutral kaon system, including the asymmetric coupling between 
two oscillation modes which leads to the phenomenon of CP violation. This system 
is the Foucault pendulum. In the present article we explore the parallels between the 
Foucault pendulum and the neutral kaon system, showing how one might construct 
a pendulum with the closest possible relation to the actual short-lived and long-lived 
neutral kaon states, and discussing not only CP, but T and CPT violation as well. 

1 To be submitted to Am. J. Phys. 
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This article is arranged as follows. We first recall the motion of the pendulum in 
Section II, review the neutral kaons in Section III (partly a recapitulation of results 
from Ref. 0), and explore the parallels in Section IV. Some suggestions for modifying 
the basic pendulum to make it emulate the actual kaon system are made in Section 
V, while Section VI concludes. An Appendix deals with the CPT-violating case. 

II. FOUCAULT PENDULUM 



A. Equations of motion in a rotating frame 



The motion of a body subject to a force F in a system rotating with constant 
angular velocity O e is described by [7] 

mr = F — 2m(f2 E x r) — mO E x (fi E x r) — 2mf3r , (1) 

where the second term is the Coriolis force, the third term is the centripetal acceler- 
ation, and the last term has been inserted to describe damping due to air resistance. 
For the case in question, O e will be a vector pointing toward the Earth's North Pole, 
with magnitude 2n d -1 . 

Consider a pendulum with a spherically symmetric support point so that it is free 
to move in two directions. Denote the corresponding axes by x (East), y (North), 
and z (up, i.e., perpendicular to the surface of the Earth). The components of f2 E 
are 







f2 E cos 9 , f2 Ez = f2 E sin 9 



(2) 



where 9 is the latitude (positive for North latitude). The components of O e x r are 
(O e x i) x = — fi E sin 9y , (f2 E x r) y = fi E sin 9x , (O e x r) z = — fi E cos^x , (3) 



where we have neglected z for small oscillations of the pendulum. 

The centripetal acceleration term — mf2 E x (f2 E x r) in Eq. (|l|) has magnitude 
mf2 E 2 a cos 9 and acts in the direction — i:cos# + ysin#, where a is the radius of the 
Earth and f2 E 2 a = 3.38 cm s~ 2 0. It thus changes the local acceleration of gravity 
slightly, g — > g c s, and displaces the equilibrium position of the pendulum. We shall 
take account of these effects by redefining g = g c g and setting x = y = to be the 
equilibrium position. 

We then define Uq = g/l, where I is the length of the pendulum, and Q = f2 E sin 9, 
write Eq. ([!]) in component form, and cancel a factor of m. The result is 



x 



—Jix + 2Qy — [3x 



y 



ujly — 2Qx — j3y 



(4) 



The coupled equations (f|) can be solved for periodic motion by substituting x = 
xq exp(— iut), y = yoexp(-iut), leading to an eigenvalue equation for u 2 . Expanding 
around uj = u and dividing by 2u , the result is 



M 



x 




x 


= UJ 









(5) 
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where 



M 



ujq — i/3 —ifl 
iQ ujq — if3 



(6) 



Eq. (H) is very close to the result which one would obtain for mixing of a neutral-kaon 
system in which the K° is represented by the basis vector x while the K° is represented 
by the basis vector y. We shall explore this parallel presently. First, however, we show 
that Eq. (|B|) leads to the familiar behavior of the Foucault pendulum in which the 
plane of linear oscillations precesses by a daily amount which depends on the latitude. 



B. Solution of equations of motion 

The normalized eigenmodes of the system (H) are 



(7) 



with eigenvalues 



fi R = uj - i(3 + Vt , n L = u -if3 -VL 



An arbitrary two-component solution x(t) describing motion in the x — y plane can 
then be written as a linear combination of \R) and \L) as 



x(t) =c R \R)e-^ t + c L \L)e 



ifi L t 



(9) 



The initial conditions on x(0) and x(0) then permit us to specify the complex quan- 
tities cr and cl, by imposing the condition 



Re x(0) 



x(0) 
2/(0) 



Re x(0) 



x(Q) 

y(o) 



(10) 



on the solution @ at t — 0. 

Let us now assume the pendulum is initially displaced in the x (East) direction, 
with zero velocity, so that x(0) = x , x(0) = y(0) = y(0) = 0. We then find 



Re c R 



x u>o — 
Im cr = Im cl 



Re cl 



xq ujo + 
V2 



(3x Q 



\/2a; 

The solution for the motion of the pendulum as a function of time is 



(11) 



'12) 



x(t) = x e pt 
y(t) = xoe-* 



(cos coot + — sin uot) cos Qt H sin Lurf sin Vtt 

(cos uj t + — sin u t) sin ilt sin u t cos ilt 

UJ Uq 



(13) 
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For fi <C u , one sees by making the replacement Sit — > fit + | and comparing the 
terms ~ sin fit and cos fit in Eq. ( |13D that the plane of oscillation rotates with angular 
frequency — fi = — fig sin 6. For = there is no precession of the plane since at 
the Equator the only effect of the Earth's rotation to lowest order is a change in the 
effective value of g. At the North or South pole the pendulum oscillates in a fixed 
plane in an inertial frame, and thus its plane of oscillation rotates with respect to the 
frame fixed with respect to the Earth with angular velocity =FfiE- 

III. THE TWO-STATE KAON SYSTEM 

The neutral kaon K° was first identified in cosmic radiation in the late 1940's ||, 
via its decay to a pair of charged pions tt + tt~ in a cloud chamber. It was "strange" 
because its production occurred much more rapidly than its decay. Gell-Mann and 
Nishijima [[J explained this feature by assigning the K° a "strangeness" quantum 
number 5=1, conserved in production processes but violated in decays. A typical 
production process, for example, would be 

7T- (S = 0)+p (5 = 0) -> K°(S= 1)+A{S = -1) , (14) 

where A is another "strange" particle first observed in the late 1940's. Strong inter- 
actions were assumed to conserve strangeness, while weak interactions (such as those 
responsible for K° — > tttt) could violate it by a maximum of one unit. 

The strangeness scheme implied that the K° could not be its own antiparticle (in 
contrast to some other neutral particles like the photon and the neutral pion tt°) , since 
it carried an additive quantum number 5 = 1. There then had to exist a neutral kaon 
with S = —1, the K . When Gell-Mann proposed this scheme, Enrico Fermi asked 
him what would distinguish the K° from the K , since they would have equal masses 
and each could decay to tttt. The solution, proposed by Gell-Mann and Pais, |§, was 
that one linear combination of the K° and K , namely K\ = (K° + K ) / a/2, would be 
able to decay to tttt, while the orthogonal combination, K 2 = (K° — K°)/\^2, would 
be forbidden from decaying to tttt and thus would be longer-lived. This particle was 
searched for and found |10|| . Its main decay modes are 3tt, ireis, and ~k\w. These 
three-body decays have smaller phase space and thus occur with a smaller rate than 
the 2tt process. 

The Gell-Mann - Pais scheme was originally based on the assumed invariance of 
the weak (decay-causing) interactions under parity (P), or mirror reflection. When 
the weak interactions were found in 1957 to violate parity invariance, the argument 
was recast in terms of the product CP of charge reflection (C) and parity, since CP 
was thought at the time to be the symmetry obeyed by the weak interactions. The 
K° and K have spin zero. A spin- zero final state of tttt necessarily has CP eigenvalue 
equal to +1, since in the center-of-mass frame inverting space and reversing the signs 
of all charges restores the original system. Thus, if CP is conserved, it is the CP-even 
linear combination of K° and K which decays to tttt. With a phase convention for the 
neutral kaons and the charge-conjugation operator chosen such that CP\K°) = \K ), 
this is just the combination Ki defined above. 
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Thus, in the limit of CP conservation, the neutral kaon states with definite mass 
and lifetime are 

* 1= *^ , K 2 = ^K . (15) 

These may be emulated with a pair of pendula of equal frequencies by coupling them 
with a spring which also dissipates energy. One pendulum may be thought of as the 
K° and the other as the ~K° . In the absence of coupling, there are two degenerate 
modes for which any two orthogonal linear combinations of K° and K can serve as 
an acceptable basis. When the pendula are coupled to one another via the spring, 
however, the linear combination in which the pendula are moving with equal and 
opposite displacements will be shifted in frequency and damped. This will be the 
mode Ki defined above if the K° and K are defined with appropriate phases. The 
mode in which the pendula are moving in phase with one another will correspond 
to the K 2 and will be unshifted in frequency and undamped. The common damping 
term for K° and K corresponding to air resistance will imply that the mode with 
the pendula moving in phase with one another is also damped, but much longer-lived 
than the mode which excites the spring. Thus, damping of both eigenmodes permits a 
meaningful parallel with the neutral kaon system. A simple laboratory demonstration 
H can be constructed which illustrates the problem up to this point. 

The coupling and damping of a two-state system in quantum mechanics can be 



described by a "mass matrix" M. very similar to that discussed in Sec. II A fllT 

: M 



d 



K° 
K° 



K° 
K° 



M = M- iV/2 , (16) 



where M and V are Hermitian. The eigenstates of M. evolve in proper time t as e~* M4 *, 
where \i{ (i = 1,2) are the (complex) eigenvalues of M.. Each /ij may be expressed 
as /ij = rrii — iji/2, where m, is the mass of the eigenstate % and 7* denotes its decay 
width. Here and in what follows we are using units with H — c — 1. The lifetime of 
the eigenstate i is then n = 1/7*. 

In the CP-invariant case when the eigenstates of A4 are given by (|l^), a little 
algebra shows that one must have 

M U = M 2 2 ; M 2 i=M 12 . (17) 

This form is only consistent with the case of the Foucault pendulum discussed in 
Sec. II when the Q term in Eq. (0) vanishes, i.e., when there is no rotation. We will 
see that rotation of the reference frame in classical mechanics is analogous to CP 
violation in the neutral kaon system. 

When CP is violated 0, the states of definite mass and lifetime K$ (for "short- 
lived") and Kl (for "long-lived") are eigenstates of a more general mass matrix no 
longer obeying A4 2 i = M.\ 2 . However, it turns out that invariance under the product 
CPT, where T denotes time reversal, still guarantees M.\\ = M. 22 - In words, this 
just says that the quantum-mechanical amplitudes for K° — > K° and K° — > K° are 



equal. CPT invariance is a very general property of quantum field theories [12fl and 



5 



we shall assume its validity here for the moment. We shall see in the Appendix how 
to relax the constraint = M. 22 . 

With CPT invariance assumed but CP violated, the eigenstates of the mass matrix 
M. can then be written as 



\S) 



2(1 + N 2 ) 
1 



2(1 



(1 + e)\K°) + (1 - e)\K°) 



{l+e)\K°)-(l-e)\K°) 



(18) 
(19) 



where S, L denote K s , K L , and e is related to M. by 



VM^ - VM^ M 12 - M 2 i Im(ri 2 /2) + i ImM; 



.12 



VM 12 + VM 2 i A^/M l2 M 2 i Us - 

Here the eigenvalues [is,L of M. are related to its elements by 



(20) 



l^s = Mn + \J M 12 M 21 ; il l = Mn - V M n M 2 i ■ 
The eigenstates \S) and \L) can be parametrized approximately as 
\S)~\Ki) + e\K 2 ) , \L)~\K 2 )+e\Ki) . 



:2ii 



:22i 



These states have a scalar product (L\S) ~ 2 Re e, which is no longer necessarily 
zero. 

Both the magnitude and phase of e can be measured accurately. The squared 
magnitude of e is measured by comparing the rates for Kl — > tctt and Ks —* tttt; the 
result is |e| ~ (2.28 ± 0.02) x 10 -3 . The phase of e can be measured by experiments 
in which decays of Ks and K^ produced in a known relative phase add coherently; 
the result is Arg e ~ 43° ||13| . Given what we know about the masses and lifetimes 

lol) 



of Ks and Kl [[13], this phase agrees with what one can derive from Eq. (|20D |lj 
Specifically, one can show \TA\ that |Imri2/2| -C |ImMi2|. This result then implies 
that 

Aige^\ 27q0 \ - Arg(/i S - [II) for <| ™ n / 2 > n \ (23) 



The observed masses m s ^ and decay widths js,L i n f^s,L = m s,L — ^7s,l/2 are such 
as to imply j§ Arg e = (43.5 ± 0.1)° for Im M n < 0, and Arg e = vr +'(43.5 ± 0.1)° 
for Im M12 > 0. 

Equation fl20| ) implies that e 7^ arises from A^i2 7^ A^2i- As we have seen, the 
Foucault pendulum with Q 7^ in Eq. (§) provides a mechanical illustration of this 
behavior. In the next section we shall explore this parallel a bit further. First, we 
discuss one quantity sensitive to e in the neutral kaon system. This is the asymmetry 
in the semileptonic decay of a neutral kaon beam to those final states which can arise 
from K° or K°. We introduce this topic because it will turn out to have a parallel in 
the case of the Foucault pendulum. 
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We shall assume, in accord with predictions at the quark level, that the allowed 
processes are K° — > ix^t+vi and K° — > 7r + £~z/^, where £ = e or //. Suppose, for 
example, that all the Kg in a neutral kaon beam have decayed away, leaving pure 
Kl- The leptonic asymmetry 

is then just ylf = 2 Re e. Measurements of this quantity are consistent with the 
measured magnitude and phase of e. 

IV. PARALLELS: BASIC PENDULUM 

The eigenstates R and L of the Foucault pendulum problem, which are eigen- 
vectors of fl6|), can be put into correspondence with neutral-kaon eigenstates S and 
L by the correspondence R *->■ S, L L. Aside from overall phases assigned to 
the states, one then sees from Eq. (p0|) that e = — i sgn(fi). The eigenstates have 
different oscillation frequencies, /is — jiL — 2fi, but their lifetimes are the same. Since 
e is imaginary, the eigenstates S and L are orthogonal to one another, just as in the 
CP- conserving case. 

The 7T7T mode is excited by the decay of K\ in (pT5|), which corresponds in the 
Foucault pendulum case to an oscillation along the line x = y. Both the eigenmodes 
(0) have a component along this direction. Since the plane of oscillation of the 
Foucault pendulum is always rotating (as long as Q ^ 0), the excitation of the irn 
mode will undergo variations in time as this plane rotates with angular frequency fl. 

The analogue of the leptonic asymmetry At discussed at the end of the previous 
Section is the asymmetry of the projection of oscillations of the L eigenmode onto 
the x K° and y <-> K° directions. As one sees from (0), there is no asymmetry, 
only a phase difference, with respect to oscillations in the x and y directions. The 
same conclusion can be drawn from the fact that Re e = in this example. 

V. MODIFICATIONS TO ILLUSTRATE ACTUAL KAON SYSTEM 

One needs eigenmodes with vastly different lifetimes in order to emulate the true 
neutral-kaon system, since T S /T L = t l /t s — 579. It is not difficult to simulate 
such eigenmodes in a two-state system [|| [|. For example, as mentioned in Sec. Ill, 
two pendula of the same natural frequency can be coupled to one another through a 
dissipative spring ||, leading to a difference in both frequency and lifetime between 
the eigenmodes in which the pendula oscillate in or out of phase with respect to each 
other. 

For the spherical pendulum, one can simply introduce damping in one of the 
two directions, for example by using a permanent magnet as the pendulum bob, 
and placing a flat coil with windings oriented along one direction just beneath the 
pendulum. The coil should be connected to a dissipative load (e.g., a resistor). A 
coil with its windings along the x direction will damp oscillations in the y direction. 
However, if we were to treat the x and y damping differently, we would be violating 
CPT invariance, since then M22 -Mil- 



7 



In the present Section we wish to emulate a CPT-preserving system, which is most 
convenient in a new basis. In the (K°,K°) basis a CPT-preserving mass matrix has 
the form 

" Mu Mvi 
M 2 i Mu 



M 



(25) 



Thus one must have equal natural frequencies and damping terms for oscillations 
in the x and y directions. However, one can transform [plj to the basis (Ki,K 2 ) 
corresponding to oscillations in the (K° ± K°) / y/2 = (x ± y)/ \[2 directions. In this 
basis the mass matrix is Af = UMW, where 



One finds 



V2 



Af 



1 



(26) 



Mi N12 
-M12 A22 



(27) 



where jV n = M n + (M u + M 2 i)/2, N 22 = Mu - (M 12 + M 21 )/2, and M2 
(M21 — Mi 2 )/2. Thus, in this basis, one can have a mass matrix of the form 



Af 



Mi — Z7x 
—ifl' uj 2 — 



(2f 



where Q' is not necessarily real. The CPT-invariance is expressed in the K\-K 2 basis 
by the condition M 2 \ = —N\ 2 . 

Physically one can emulate the system described by the matrix M by having a 
Foucault pendulum with different damping in the K\ and K 2 directions. Such damp- 
ing could be implemented, for example, by the inductive setup noted above. Natural 
frequencies in two orthogonal directions can be made to differ using a hinged set of 
supports. One could also introduce different damping constants in two perpendicular 
directions through properties of the hinged joints themselves. 

The Foucault pendulum case corresponds to real Q' for the new basis. Note that 
Q' in Af is then the same as Q in M; aside from a sign flip in the off-diagonal elements, 
Af and M are the same matrix. 

For loi — iji ^ uj 2 — 272 and fl <^ uji i2) the eigenvectors and their corresponding 
eigenvalues are approximately 



\S) 



1 



-in 



Us = loi- Z7 X + 8 S , 

n 2 



LO1- U2- 2(71 - 12) 



Ss 



U1- U2 - 2(71 - 72) 



\L) 



1 



-m 



LO1-L02- 2(71 - 72) 



\i L = L02 - 112 + 8 L , 
-- e s = e , 5 L = -ine = -5 S 



(29) 

(30) 
(31) 
(32) 
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Let us now discuss the time-evolution of states which are initially \S) and \L): 



\S) 



' 1 " 




' 1 " 


e 




e 



-ifist 



\L) 



(33) 



Defining = Arg e, these lead to the time-dependences 
Re xs(£) = x Re 



cosa;it 
e| cos(u;i£ — 0) 



Re xx,(t) = xoRe 



-72* 



|e| cos^t - 
cos a^i 



(34) 



(35) 



where contributions of order Ss,l have been omitted. Thus both S and L eigenstates 
correspond to orbits which are slightly rotated to favor an orientation in the direction 
of K°-\ike decays by an amount proportional to e. The Coriolis force in this example 
induces an effect which is akin to the leptonic asymmetry parameter Ae discussed 
at the end of Sec. III. Moreover, the decay to tttt in the K\-K 2 basis is like the x- 
projection of the eigenstate. Although the dominant damping is along the x-direction 
(we assume 71 72), the presence of the Coriolis force causes the eigenstate L to 
have a projection proportional to e along x, i.e., the K L does decay to nn. 



VI. CONCLUSIONS 



We have shown that the motion of a Foucault pendulum has many features in 
common with the CP- violating neutral kaon system. When the natural frequencies for 
oscillation in two perpendicular directions and the damping terms for these directions 
are equal, the parameter e describing the eigenstates \K S ) — \Ki) +e\K 2 ) and \K L ) ~ 
\K 2 ) + e\Ki) takes on the special value e = —%. In order to simulate the physical 
situation in which |e| = 0{2 x 1CT 3 ) and Arg(e) ~ 7r/4, one must construct a Foucault 
pendulum with slightly different natural frequencies in two perpendicular directions, 
and with vastly different damping constants in these directions. While the practical 
realization of such a construction sounds challenging, it is interesting that, at least in 
principle, it seems feasible entirely within the realm of classical mechanics. 

As we show in the Appendix, the phenomenon of CPT violation, with preservation 
of T and violation of CP, can be emulated by coupled resonant circuits, building upon 
the results of Ref. Q. The program set forth in that work is still incomplete; we would 
be delighted to see an implementation of CP and T violation, with CPT conservation, 
through the asymmetric coupling of two resonant circuits with equal frequencies. 

The classical emulation of CP violation via the Foucault pendulum leaves us with 
one big puzzle. In the classical system, the asymmetry in coupling between the two 
modes is imposed from the outside, so to speak, through the Earth's rotation. In the 
neutral-kaon system, the corresponding asymmetry in K°-K° mixing is thought to 
arise from a complex phase in the Cabibbo-Kobayashi-Maskawa matrix |TJ| describing 
the charge-changing weak transitions of quarks. Is that phase an indication of a new 
fundamental asymmetry arising from physics beyond the Standard Model? 
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APPENDIX. CPT- VIOLATING CASE 

The matrices M. and M are arbitrary when CPT is violated. The relation between 
them is 

Mi = {Mii+Mii + Mn + M&)/2 , Afi2 = {Mii-Mi2 + M2i-M 2 2)/2 , 

N21 = (M 11 +M 1 2-M 2 i-M 2 2)/2 , JV 2 2 = (Mn-Mi2-Mn+Mvi)/2 .(36) 

A simple example of a CPT-violating mass matrix M involves coupling between two 
resonant circuits, as discussed in Ref. j|. If the circuits are taken to have different 
frequencies, the matrix takes the form 



Af 



uj\ — 271 a 
a uj 2 - h2 



(37) 



in which the off-diagonal elements are equal (rather than equal and opposite as in 
the CPT-preserving case). An analysis parallel to that for the eigenstates S and L 
performed in the previous Section leads to the results (for \a\ <C 00%^) 



\S) 



1 



1 



with 



a 



U1-UJ2- - 72) 



-e L = e, 



a 



fi S = ui-i f yi-\ 

LU1-U2- H7i - 72) 

The eigenstates are thus 



, /i L = ^2-272 • 



a 



LU X - LU 2 - £(71 - 72) 



|S)~|l) + e|2) , |L)~|2)-e|l) 



(38) 

(39) 
(40) 

(41) 



This case (see, e.g., Ref. |[16|| ) correponds to invariance with respect to time-reversal, 
so that CP and CPT are violated. Both ( p2]) and ( T4T| ) can be written in the more 
general form 

|5)~|l) + e s |2) , |L)~|2) + e L |l) . (42) 

The CPT-preserving, CP- violating case corresponds to es = = e, while the case 
fl41"D corresponds to = — e L = e. To lowest order in es l, one finds 



\K ) = -L[\S)(l-e L ) + \L)(l-e s )} , \K°) = -L[\S)(l + e L )-\L)(l + e s )} .(43) 
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Since the states \S,L) evolve with proper time t as \S, L) — > e lfls > Lt \S, L) , a short 
calculation shows that 

\K°) -> \K°)[f + (t) + (e s - e L )/_(t)] + \K°)[1 - ( e5 + e L ))f4t) , 

\K°) -> + (<fc - + |A-°)[1 + (6 5 + e L )]/_(t) , (44) 

where f±(t) = [exp(— ifist) ± exp(— i/iLt))/2. 

For = —65, the evolution of K° into A' is the same as that for K° into K°, 
corresponding to a time-reversal-invariant situation. However, the amplitudes for 
K° — > K° and A' — > A' differ from one another, corresponding to CPT violation. 

For 6l = 65 = e, the amplitudes for K° — > X° and K° — > K° are the same, 
corresponding to CPT invariance, while those for K° K° and K° — >■ X° differ 
from one another, corresponding to T violation. In this case the terms |/_(t)| 2 cancel 
in the rate asymmetry 

A = r[^(o)^Z (t)]-r[Z°(o)^K (t)] 

T T[K°(0) ^K°(t)} + T[K°(0) AT°(t)] ' 

and to lowest order one finds |T^] ^4t = 4 Re e. This relation has recently passed an 
experimental test at CPLEAR pR . 
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